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ABSTRACT 

IC 4051 is a giant E galaxy on the outskirts of the Coma cluster core. 
Using archival HST WFPC2 data, we derive the metallicity distribution, 
luminosity function, and spatial structure of its globular cluster system (GCS). 
The metallicity distribution derived from the (V — I) colors has a mean 
(Fe/H) ~ —0.3, a near-complete lack of metal-poor clusters, and only a small 
metallicity gradient with radius. We tentatively suggest that the GCS has 
two roughly equal metallicity subcomponents, one centered at [Fe/H] ~ 0.0 
and the second at [Fe/H] ~ —1.0, although their identification is blurred by 
the photometric uncertainties. The luminosity distribution (GCLF) has the 
standard Gaussian-like form observed in all other giant E galaxies, with a 
peak (turnover) at V° = 27.8, consistent with a Coma distance of 100 Mpc. 
The radial profiles of both the GCS and the halo light show an unusually 
steep falloff which may indicate that the halo of this galaxy has been tidally 
truncated. Lastly, the specific frequency of the GCS is remarkably large: we 
find Sn = 11 ± 2, at a level which rivals M87 and most others in the central 
cD-type category, even though IC 4051 is not a cD or brightest cluster elliptical. 
This galaxy exhibits a combination of GCS characteristics found nowhere else. 
A formation model consistent with most of the observations would be that 
this galaxy was subjected to removal of a large fraction of its protogalactic gas 
shortly after its main phase of globular cluster formation, probably by its first 
passage through the Coma core. Since then, no significant additions due to 
accretions or mergers have taken place, in strong contrast to the central Coma 
galaxy NGC 4874. 



Subject headings: Galaxies: Star Clusters; Structure; Evolution; Elliptical; 
Invidual 
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1. INTRODUCTION 

The Coma Cluster (d ~ 100 Mpc), as a rich Abell cluster, is the host environment for 
a huge range of E/SO galaxies. Its two central supergiants, NGC 4874 and 4889, are among 
the very most luminous galaxies known, and the cluster has many other large ellipticals 
scattered throughout its ~ 1— Mpc core region and well beyond. The globular cluster 
systems (GCSs) around these galaxies are well within reach of the HST cameras, and thus 
give us an extraordinarily rich range of target galaxies for comparative GCS studies. Coma 
also has important implications for the determination of the Hubble constant Hq, since 
its large recession velocity (~ 7100 km s -1 ) greatly exceeds any anticipated local peculiar 
motions, and it is situated at high galactic latitude {b = +87? 7) nearly unaffected by 
foreground absorption, 

IC 4051 is a giant E2 galaxy located ~ 14' east in projection from the center of Coma 
and has no luminous neighbors. With an integrated magnitude = 13.20, it is the fifth 
brightest elliptical galaxy within the 1° central region of Coma. Thus, IC 4051 presents a 
good opportunity to study the GCS of a more-or-less "normal" large elliptical. Paum et~aT 



( 1997)| obtained the first deep WFPC2 exposures of IC 4051 with exactly this purpose in 



mind, but they limited their published analysis to only the photometry from the PCI frame 
on which the galaxy was centered. However, considerably more information on the radial 
structure of the GCS, its specific frequency, and radial metallicity gradient are potentially 
available from the outer sections of the WFPC2 field. In this paper, we re-analyze the 
archival HST/WFPC2 frames from the Paum et al. (1997]1 program and discuss the global 



properties of the IC 4051 GCS using the entire body of data available in the raw exposures. 



2. DATA ANALYSIS 

The raw database comprises HST/WFPC2 observations taken on 1995 July 20/21 
(Baum et al. 1997). Eight images totalling 20500 s (essentially, a sequence of full-orbit 
exposures) were taken through the FQ06W (wide V) filter, and two images totalling 5200 s 
were taken through the F814W (I) filter. The techniques in our data analysis are similar 
to those in Papers I and II of this series ( |Kavelaars et al. 2000| ; |Harris et al. 2000| ), in 



which we analyzed the GCS around NGC 4874, the central cD giant. The first step was 
to register and combine the images to produce a master V and master / frame free of 
cosmic rays. An excellent color composite image constructed from the combined frames is 
published by Baum et al. For the PCI frame, we produced a "flattened" image in which the 
overall elliptical contours of the central galaxy were modelled and then subtracted. In the 
WF2,3,4 fields, we generated an empirical model of the galaxy light by median filtering each 
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frame and subtracting the smoothed image from the original picture (after a preliminary 
star- finding and removal). 

Our photometry was then carried out on the flattened master frames with the 
DAOPHOT II and ALLSTAR codes ( (Stetson 1992|) . The instrumental magnitudes returned 



by ALLSTAR were transformed to the standard V,I system via the equations derived by 
|Holtzman et al. (1995) 



V = m(F606W) - Am ± 22.093 + 0.254(V - I) + 0.012(V - I) + 2.5 log(GRi), (1) 

I = m(F8UW) - 5m ± 20.839 - 0.062(^ - I) + 0.025(y - I) 2 + 2.5 log(G^), (2) 

where m(F606W, F81AW) are the instrumental magnitudes returned by ALLSTAR, and 
the Am are constants which shifts these ALLSTAR psf magnitudes to the equivalent 
magnitude of the light within a 0"5 aperture. The mean value of Am was determined 
empirically for each of the four WFPC2 fields from roughly ten moderately bright, isolated 
stars in each CCD, and applied to all our detected starlike images. The typical internal 
uncertainties in the mean Am values were ±0.02 mag. 

It should be noted that at the distance of Coma, all globular clusters are easily starlike 
in appearance and thus these transformation equations normalized to the 0"5 aperture can 
strictly be applied to them. The calibration equation for V quoted above is the synthetic 
model transformation from Table 10 of |Holtzman et al. (1995)| , while the equation for I is 



the observationally based one from their Table 7. If, instead, we had chosen their Table 10 
model transformation for J, our resulting (V — I) color indices would have ended up bluer 
by roughly 0.02 mag at the average color (V — I) ~ 1.2 of the IC 4051 globular clusters 
(see below; this shift would have the effect of reducing the deduced cluster metallicities by 
A[Fe/n] = 0.1).[] When we add the internal uncertainties in the aperture corrections Am, 
we estimate that the zeropoints of either the magnitude or color scales are uncertain to at 
least ±0.03 mag. 

Finally, as an internal check on the relative zeropoints of the V and I scales between 
the PCI chip and the three WF chips, we inspected the mean (V — I) color indices of the 
measured globular clusters in the annulus around the center of IC 4051 (R ~ 10" — 20") 
that overlapped the PC/WF boundary. We enforced the WF2,3,4 data in this annulus to 



1 Baum ct al. (1997) adopted the synthetic transformation for / from Table 10 of Holtzman et al., but 
used the specific curve for the color range (V — I) < 1.0. However, the majority of the globular clusters are 
somewhat redder than this, so that the transformation for (V — I) > 1.0 would have been preferable. The 
net result is to increase the small (V — I) difference between their scale and ours by a further ~ 0.02 — 0.03 
mag, although our instrumental color scale (F606W — F814W) agrees quite closely with theirs. 
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have the same mean color as those in the PCI region by adjusting the I magnitudes (which 
are from much shorter exposures than V, thus internally more uncertain). The final result 
places our (V — I) scale in agreement with the [Baum et al. (1997}] scale to within 0.03 mag, 
a level entirely consistent with the combined photometric uncertainties. 



As was already evident from the Paum et al. (19971 study, IC 4051 has a populous 



GCS, which appears as a very obvious swarm of faint objects across the whole WFPC2 field. 
The main source of sample contamination is from very faint, compact background galaxies, 
with a (nearly negligible) contribution from Galactic foreground stars. A high proportion 
of the background galaxies can be eliminated through conventional radial-moment image 
analysis. For this purpose we used the r\ radial moment as implemented in [Harris et al 



n = , (3) 

which is an intensity- weighted mean radius for the object calculated over all pixels brighter 
than the detection threshold (see Harris et al. 1991). A straightforward plot of r\ against 
magnitude then shows a well defined stellar sequence, with nonstellar objects scattering 
to larger r±. These classification graphs for the four CCDs in V are shown in Figure 
The dashed lines indicate the adopted cutoffs applied to the measurements. Similar object 
classifications were applied to the / data (which, however, have shallower limits), with the 
final culled data lists containing 4058 objects in V and 1672 objects in (V — I). 

The faint-end completeness of our photometry was investigated through an extensive 
series of artificial-star tests on the master images. The procedure performed here was to add 
500 artificial stars to each frame over a range of input magnitudes, measure these frames 
through the normal DAOPHOT sequence, and find out how many were recovered. Fifteen 
of these trials were carried out, with average resulting completeness fractions as displayed 
in Figure 0. Convenient fits to the raw points are provided by the Pritchet interpolation 
function, 

a(R — Rum) 



l + a 2 (R-R, 



lim I 



\2 



(4) 



Table [I] summarizes the best-fit parameters to Equation [| for each CCD and bandpass. In 
the Table, (Vu m , Iu m ) are the magnitudes at which / drops to 0.5, and the parameter a 
controls the steepness of the falloff. The curves for all three of the outer chips (WF2,3,4) 
are nearly identical; for the inner PCI chip, the limiting magnitudes are brighter, driven by 
the spread of the PSF over many more pixels and (for R< 10") the brighter background 
light. 



Lastly, Figure [| shows how the photometric measurement uncertainties (also derived 
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from the ADDSTAR completeness tests) increase with magnitude. At the formal limiting 
(/ = 0.5) magnitude, the rms uncertainty in the photometry reaches 0.15 mag. Wherever 
possible, we avoid dealing with any features of the data below that limit. 



3. COLOR AND METALLICITY DISTRIBUTIONS 

The distribution in colour of the globular clusters can be used to gain insight into the 
existence of multiple sub-populations in the GCS. Bimodal color distributions are found 
about half the time in gE galaxies (e.g., |Kundu fc Whitmore 199§| pNeilsen fc Tsvetanov 



1999[ ) and are often interpreted as relics of at least two major phases of star formation in 



the early history of the galaxy, whether by merger, accretion, or in situ processes. With 
the conventional "null hypothesis" for giant ellipticals that the clusters are all old ( > 10 
Gy), the color index is primarily a tracer of cluster metallicity. The colour- magnitude 
distribution for the 1672 objects measured in both V and I is shown in Figure f|. At 
projected galactocentric radii larger than about 80", we found (see below) that the residual 
numbers of clusters dropped nearly to zero, so we adopt this outer region as defining a 
suitable "background" population. 

To eliminate a few more contaminating objects, we reject objects bluer than 
(V — I) = 0.74 or redder than (V — I) = 1.46 (vertical dashed lines in Fig. |j). These 
colour limits generously include the range in colours of the known globular clusters in 
large galaxies (e.g., [Harris 1996| ; [Whitmore et al. 1995| ; |Neilsen fc Tsvetanov 19991) . We 



also further limited the color sample to objects brighter than V = 26.0 to ensure high 
completeness at all colors. 

A "clean" color distribution for the GCS was then obtained by subtraction of the 
background (R > 80") color distribution, normalized to the same total the inner 

population. This procedure left a final total of 479 objects within the magnitude and color 
limits given above, with a net distribution over (V — /) as shown in the histogram of Figure 
IH The mean color of the sample is (V — I) = 1.12 ± 0.01 (internal uncertainty of the 
mean), with a dispersion of oy_j = 0.13. Subtracting an adopted foreground reddening 
E(V — I) = 0.014 and using the calibration of (V — I)o in terms of metallicity given in 
Paper II, 

(V-I)q = 0.17[Fe/H] + 1.15, 

we then estimate that the IC 4051 GCS as a whole has (Fe/H) ~ —0.3. The peak position 
of this color distribution is quite similar to the metal-rich components in other giant 
ellipticals such as NGC 4472 QGeisler et al. 1996| ), M87 QWhitmore et al. 1995fc |Kundu et 



|al. 1999| ), and other Virgo members ( pNeilsen fc Tsvetanov 199"9] ). However, the metal-^oor 
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component which is usually found in these same galaxies at a mean color (V — I) ~ 0.95 or 
[Fe/H] ~ — 1.5 (and which we found in the Coma cD NGC 4874; see Paper II) is entirely 
missing in IC 4051, or at very most is a fringe component buried in the wings of the main 
distribution. 

We cannot place firm limits on the intrinsic dispersion <r[Fe/H], since the mean 
observational measurement scatter over the sample is uy-i — 0.11, comparable with 
the observed sample dispersion of ±0.13. Nevertheless, subtracting off the observational 
scatter in quadrature, we estimate roughly o"o[Fe/H] ~ 0.4, which is in close agreement (for 
example) with the value cr[Fe/H] = 0.38 found by Geisler et al. (1996) to fit each of the 
metal-rich and metal-poor components in NGC 4472. In the Milky Way, the well known 
bimodal MDF has been found to be fit by Gaussian functions with dispersions near 0.3 
dex ( |Zinn 19E5] ; [Armandrotf fc Zinn 1988] ; |Harris 2000| ). For IC 4051, a single Gaussian 



with the same mean and standard deviation as the sample (Fig. ||) matches the MDF with 
a x 2 — 14.6 over 14 degrees of freedom, which provides no strong evidence for bimodality 
(but see below). 

Trends of mean color with either galactocentric distance or magnitude were also 
searched for. Table |2| shows the mean color and dispersion in 0.5-magnitude bins from 
V = 22.5 to 26.0. These binned means reveal no significant change in color with luminosity. 
However, slightly more interesting features emerge in the graph of color versus radius 
(Figure ^). Binned mean colors, listed in Table || indicate no systematic change in color 
for R> 10", but within 10" the clusters are indeed slightly redder than the overall mean. 
The distribution in its entirety is barely suggestive of two sub-populations: one centered on 
(V — I) ~ 1.2 which is found at all radii; and a second, slightly bluer one centered near 
(V — I) ~ 1.0. The lack of bluer clusters within R< 10" is then largely responsible for the 
inner color gradient of the whole sample mentioned above. Much stronger versions of this 
same effect have shown up in some other giant E or cD galaxies with far more obvious 
bimodal MDFs (e.g., [Seeker et al. 1995[ |Geisler et al. 1996| ; |Lee et al. lM| |Ustrov fc 



[Forte 19981 ). In these, the different central concentrations of the metal- rich and metal-poor 
subsystems produce a steady outward change in the relative proportions of blue-to-red 
clusters with radius and thus a mean metallicity gradient. 

Using Fig. || as a guide, we divided the sample of objects at (V — I) — 1.07 and tested 
the radial distributions of the bluer and redder halves. A standard Kolmogorov-Smirnov 
two-sample test indicated that their spatial distributions are significantly different (the 
redder half is more centrally concentrated) at the 99% level, suggesting to us that the inner 
gradient is indeed a real effect. 

We therefore very tentatively suggest that the IC 4051 system may contain a bimodal 
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MDF in which the two modes are rather closely spaced in mean metallicity, thus heavily 
blurred out by the raw photometric measurement uncertainty. Numerical experiments 
with various two-component fittings of the entire MDF lead to models of the form shown 
in Figure [F[ Here, a sample twin-Gaussian fit is shown in which the bluer (metal-poor) 
component is centered at (V — I) = 1.00 or [Fe/H] ~ —0.96, the redder (metal-rich) one at 
(V — I) = 1.17 or [Fe/H] ~ +0.04, both have dispersions cr(V — I) = 0.10, and the redder 
one contains about 55% of the total sample. The combined components now represent the 
total shape of the MDF better, with its modest skewness toward the red side (the total x 2 
is 12.6). The relative proportions of blue and red components, however, are quite uncertain 
(the formal uncertainties are ±0.1, but variations of factors of two in the proportions give 
scarcely different overall fits). 

Clearly, this particular two-component model is only illustrative of the range of 
possibilities: the moderately small difference in color between the two components, and the 
very significant broadening of the MDF by photometric measurement uncertainty, do not 
justify more extensive analysis. However, it would clearly be of value to measure the MDF 
of this populous globular cluster system with a photometric index much more sensitive to 
metallicity than (V — I), in which the subpopulations would be far more clearly revealed. 
A more sensitive color index would also permit establishment of the true mean [Fe/H] with 
much less zeropoint uncertainty. 

Lastly, it is worth comparing the mean colors of the GCS components to that of the 
halo light of the central galaxy. |Mehlert et al. (1998)| find (V — I) ~ 1.30 at a projected 



radius R = 10", increasing inward to (V — I) = 1.35 at the very center. This color range 



is distinctly redder than the typical levels (V — I) ~ 1.20 ± 0.03 for giant E galaxies ( |Buta 



|fc Williams 1995| ; [Prugniel fc Heraudeau 1998| ). The measured absorption line indices 



(Mg, Fe, H/5) lead [Mehlert et al. (1998)| to conclude, in line with the integrated color, 



that the core of IC 4051 is extremely old and very metal-rich, perhaps as high as [Fe/H] 
= +0.25. However, the deduced metallicity from the line indices becomes lower at larger 
radii, dropping to an equivalent [Fe/H] ~ —0.5 for R>20" (the effective radius r e of the 
light profile), similar to the inner GCS. Mehlert et al. find that IC 4051 harbors an old, 
co-rotating core with an unusually large "break radius" (it is detectable out to 5" or 3.4 
kpc) but which contributes < 1% of the total light of the galaxy. If this inner stellar disk is 
a signature of a dissipational merging event, it is likely to have occurred at early times. 
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4. THE LUMINOSITY DISTRIBUTION 



As |Baum et al. ( 1997}] showed, the V photometry reaches faint enough to reveal the 



"turnover point" (peak frequency) in the globular cluster luminosity function (GCLF). By 
adding in the photometry from the WF chips, we have been able to double the total sample 
of clusters and thus improve the definition of the GCLF. 

The distribution of all the detected objects classified as "starlike" and used to define 
the GCLF is shown in Figure [8[ These are, quite evidently, strongly concentrated to the 
center of IC 4051 (much more so than in the GCS of the Coma supergiant NGC 4874; see 
Papers I and II). More or less arbitrarily, we take the region R > 80" marked by the outer 
dashed line in Fig. ^ as defining the luminosity function of the background population, 
to be subtracted statistically from the inner (10" < R < 80") zone after correction for 
photometric incompleteness. 

The results of this exercise for each of the four CCD chips separately are shown in 
Figure ^. Aside from the noticeably brighter completeness limit for the PCI zone, no 
significant differences in the GCLF shape or turnover from place to place are evident. (The 
GCLF peak for the PCI region shows an apparent peak fainter than V ~ 28, but this is 
fainter than the 50% completeness limit and so cannot be given much weight.) We therefore 
add all four sectors to form the composite GCLF shown in Figure The numerical results 
in 0.3-mag bins are listed in Table £|: here, successive columns give (1) the V magnitude 
range of the bin (2) the number of detected starlike objects in the inner (10" — 80") zone 
(3) the number in the outer (> 80") background zone (4) the number in the inner zone 
corrected for completeness, and (5) the net GCLF, after subtraction of the area-normalized 
background counts. 

To estimate the turnover level and shape of the GCLF, we fit a standard Gaussian 
interpolation function ([Harris 1991| ; |Jacoby et al. 199"2|) to the data shown in Fig. 10], setting 



the standard deviation oy of the curve and then solving for the best-fit turnover level V°. 
Trials with different adopted cry's gave the results summarized in Table [|. The reduced % 2 
values favor a solution in the broad range ay ~ 1.4 — 1.8, with little to choose among values 
in this range in a formal sense. However, it is well known that both ay and V° tend to be 
overestimated in situations like these where the magnitude limit of the data reaches barely 
past the actual turnover (e.g., |Hanes fc Whittaker 1987^ Paper I) since the solutions for 



the two parameters are correlated. For this reason, we favor a choice in the narrower range 
oy — 1-4 — 1.6 and V° ~ 27.6 — 28.0. Sample Gaussian curves for the extremes of this range 
are shown in Fig. [K| 

Our final adopted pair of parameters is V° = 27.8 ±0.2, ay = 1.5 ±0.1. For comparison, 
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Baum et al. (1997) found V° = 27.72 employing a different and more complex fitting 



function. As is discussed more extensively in Paper I, this turnover level is also similar to 
what we found in the central cD NGC 4874. Using both of them combined, along with a 
calibration of the absolute magnitude of the turnover point based on the Virgo ellipticals, 
we find d ~ 100 Mpc for Coma along with a Hubble constant Hq ~ 70 (see Paper I). 

A second and more physically oriented way to display the same material is as the 
luminosity distribution function (LDF), or number of clusters per unit (linear) luminosity. 
(The relation between the GCLF and LDF forms is exhaustively discussed by [McLaughlin] 
|1994[ ) The LDF is shown in Figure O. At levels brighter than the GCLF "turnover" 
(which in turn is only slightly brighter than the photometric completeness limit), the LDF 
clearly approximates a power-law falloff toward higher luminosity, N(L)dL ~ L~ a . To 
second order, however, the slope a = — dlog(N) / dlog(L) appears to steepen slightly at the 
upper end: an unweighted least-squares fit to all bins brighter than the turnover yields 
a = 2.05, while exclusion of the half-dozen very brightest bins yields a = 1.75. 

These power-law forms - as well as logarithmic slope values a ~ 2 - are entirely similar 
to what has been found in a wide range of other galaxies from dwarf ellipticals to spirals 
( [Harris fc Pudritz 1994 ; Uurrell et al. 199T ). However, in most giant ellipticals studied to 



date, the slopes tend to be somewhat flatter at a ~ 1.5 ± 0.3 (Harris fc Pudritz 1994j ). 
The total shape for log (L/L e ) > 5, complete with its progressive steepening toward higher 
luminosity, can be well matched by a formation model in which protocluster clouds build 
up by collisional agglomeration and in which the more massive clouds have shorter lifetimes 
against star formation ([McLaughlin fc Pudritz 1996| ; [Harris 2000|) . Our data for IC 4051 



add further to the general body of material which indicates a remarkable place-to-place 
similarity in the luminosities of old globular clusters, and thus a quasi-universal formation 
process. 



5. RADIAL DISTRIBUTION AND SPECIFIC FREQUENCY 

Because the GCS around IC 4051 is quite centrally concentrated (see Fig. ||), we can 
use the complete WFPC2 data to define the spatial distribution outward nearly to its 
limits. The radial profile of the raw counts for all starlike objects brighter than V = 27.0, 
for which the data are highly complete nearly in to the central core of the galaxy, is shown 
in Figure [12|. The inner core (R<5"), in which the projected density of clusters is nearly 
flat, continues outward to a steep power-law falloff which covers most of our survey area. 
Finally, for R > 80", the number density a begins to level off towards its eventual far-field 
background level; more or less arbitrarily, we set this background at cr& = (0.02 ± 0.01) 



-10- 



arcsec -2 as representing nearly the average of the outermost two points. (As will be 
seen below, small differences in the adopted will not have major effects on any of our 
subsequent conclusions.) 

The complete profile data broken into circular annuli are listed in Table |], giving the 
number of objects in each bin, the surface area of the annulus, and the projected density 
a. The residual number density of clusters, o c \ = a — a^, is plotted in Figure [H| Simple 
King (1966)| models can be fitted to the a c i data points to give rough estimates of the GCS 



core radius and central concentration: performing a weighted fit in the manner described in 
Paper II, and ignoring the very uncertain outermost three points, we find a formal best-fit 
core radius R c = 10?25 (equivalent to 5.1 kpc at the adopted Coma distance) as well as a 
concentration index c = 1.45 for a dimensionless central potential Wq = 6.26. The core 
radius is four times smaller in IC 4051 than the ~ 22 kpc value we found in the much more 
extended NGC 4874 (Paper II). 

A second comparison can be made with the halo light of the galaxy. It is conventionally 
found in giant ellipticals that the GCS is a more spatially extended system as a whole than 
the halo ( [Harris 1991| , 1999). IC 4051 is no exception, despite its overall compact structure. 
In Fig. [H| we show the wide-field surface intensity profiles in measured by |Strom fc 



Strom (1978)| and |Jorgensen et al. ( 1992 )| . Although the Strom data are photographically 
measured, their profile agrees tolerably with the more recent CCD measurements of 
Jorgensen et al. (1992)| over their region of overlap. 



The bulk of the GCS profile is more extended than the halo light, except possibly for 
the outer (i?>30") regions where their slopes are more nearly similar. For R>20", the 
GCS profile behaves as a d ~ R~ 2 , although at the largest radii little weight can be placed 
on the very uncertain outermost half-dozen points. There is a strong hint from the halo 
light profile that the galaxy may be truncated past R ~ 60" (about 30 kpc), though here 
again the profile is very sensitive to slight differences in the adopted background level, so 
not much meaning can be ascribed to the slope differences between the GCS and the halo 
there. For giant E galaxies in general, a rough mean relation between galaxy luminosity 
My and the radial falloff outside the central core is ( |Kaisler et al. 19"96 ) <ilogcr/<ilog-R 



~ -0.29My - 8.00. For IC 4051, this relation would predict a d ~ i?" 1 ' 65 , somewhat flatter 
than the observed R~ 2 trend. 

Calculating the total GCS population and specific frequency is now a straightforward 
matter. From Table we multiply a c i(R) by the area of each annulus, then sum the annuli 
to get the total cluster population out to the limits of our survey. We find N = (1845 ± 165) 
for V < 27.0 and R< 130". If the true GCLF turnover magnitude is at V° = 27.8 ± 0.2 (see 
above), then we must multiply this raw total by ~ (3.35 ±0.52) to estimate the total cluster 
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population over all magnitudes, giving N d = 6180 ± 1100. The integrated luminosity of the 
galaxy is V T = 13.20 (RC3 catalog value), corresponding to My = —21.9 for our adopted 
Coma distance. Thus, the specific frequency is 

^ = iV c/ • 10 a4(M ^ +15) = 10.8 ±1.9. 



In strict terms this is a lower limit to the true global Sn, since we have not accounted 
for any cluster population outside the ~ 120" radial limit of our WFPC2 field. However, 
given that the halo is clearly declining quite steeply in this region (Fig. [TB|), any such 
population correction is likely to be small. A generous but reasonable upper limit estimate to 
the total population can be made if we assume that the GCS profile continues as a c i ~ R~ 2 
outward to the nominal tidal radius at Rt ~ 230". This assumption gives an additional 
~ 380 ± 300 clusters brighter than V = 27, which then translates to Sn = 12.6 ± 2.6. 

Placing more weight on the lower limit - which reflects the steep falloff of the system 
near the radial limit of our data - we adopt a final estimate 

Sjv(final) =11 ±2. 



Remarkably, this GCS population ratio is several times higher than the Sn 'Z 2 value 
found in NGC 4881 ( [Baum et al. 19951) , a galaxy which is quite comparable with IC 4051 
in luminosity, structure, and location on the outskirts of the Coma core. This high Sn, in 
fact, places IC 4051 in the range which is conventionally reserved for the central-giant cD 
galaxies like M87 and many other BCGs QHarris et al. 1998; ; [Blakeslee 1991 , 1999). It is, 
perhaps, particularly noteworthy that IC 4051 has a specific frequency three times higher 
than the central cD in its own host galaxy cluster, NGC 4874 (see Paper II). No other 
instance of such a large contrast between a low— Sn central cD and a higher— Sn outlying 
elliptical is known. IC 4051 provides striking evidence that a central location in a rich 
cluster environment is not required to form a high population of globular clusters. 



6. DISCUSSION 

A brief summary of our findings for IC 4051 is that its GCS is (a) almost entirely 
metal-rich, albeit possibly with two narrowly separated subcomponents; (b) relatively 
compact in radial structure; and (c) a "high specific frequency" system despite that fact 
that its host galaxy is not a central giant elliptical nor one with a cD-type envelope. 

Just as in Paper II for NGC 4874, we now attempt to use the integrated characteristics 
to reconstruct a partial history of the system. Formation scenarios for giant ellipticals tend 
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to fall into three basic camps: (a) "in situ" formation, whereby the galaxy condenses by 
dissipative collapse of gas clouds in its immediate vicinity, in one or more major bursts; (b) 
later mergers of pre-existing disk-type galaxies with both gas and stars; or (c) successive 
mergers or accretions of smaller gas-poor satellites. Various combinations of these extremes 
are, of course, possible, and even likely. 

For IC 4051, the lack of /ow-metallicity clusters already places fairly strong constraints 
on the range of possible formation events. For example, the mechanism investigated by 
Cote et al. (1998)| - in which an original metal-rich "seed" gE accretes dozens or hundreds 



of smaller satellites - is unlikely, since these dwarf satellites would have brought in a 
population of hundreds or even thousands of low-metallicity clusters, which we do not see. 

Similarly, merger-formation models in which gas-rich disk galaxies combine to build a 
composite elliptical QAshman fc Zepf 1992|) would predict a strong component of metal-poor 



clusters in the resulting MDF from the globular clusters that were present in the pre-merger 
galaxies. These merger models also have severe difficulty in generating high specific 
frequency products, since increasing the cluster population relative to the field stars by a 
large enough amount to produce high Sn requires very large (> 1O 1O M ) input gas masses, 
more than is routinely available in disk galaxies today. The normal merger route does 
appear to be quite effective as a logical source for low— Sn field ellipticals (see [Harris 2000] 
or [Whitmore fc Schweizer 1995] for much more extensive discussion) . 



However, if either the merger or accretion processes are taken to an extreme form in 
which the merging objects are almost completely gaseous, then they become closely similar 
to the in situ route, and the conundrum of the missing low-metallicity clusters can be more 
easily circumvented. If the gas supply - however it was assembled - underwent most or all 
of its star formation in the high-pressure, high-density environment of the protoelliptical, 
then the conversion of gas to stars would have run much further to completion and built 
up the metallicity to the high levels that we now observe. Later gaseous mergers are, of 
course, not ruled out: the central corotating disk in the core of IC 4051 (Mehlert et al. 
1998) with its very high metallicity is a likely signature of such an event, though at its < 1% 
contribution to the present-day luminosity, it probably did not form more than a few dozen 
globular clusters along with it, and even these would have mostly disrupted by now if they 
resided in the central few kpc of the core. 

The relatively compact structure of the galaxy may be the result of tidal trimming 
( "harrassment" ) from the Coma potential well (e.g., [Moore et al. 1996|) . The radial velocity 



of IC 4051 (4940 km s _1 ) is almost two standard deviations away from the Coma centroid 
(6850 km s -1 ; see |Colless fc Dunn 1996Q , indicating that this galaxy oscillates back and 
forth through the cluster and is now passing through the dense Coma core at high speed. 
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These elements of an evolutionary scenario for IC 4051 are in strong contrast to NGC 
4874, for which we argued (Paper II) that a large fraction of its clusters (which are almost 
entirely low-metallicity) could have been acquired by accretions of smaller satellites. In 
IC 4051, we are forced to argue that the bulk of its clusters formed in situ. The globular 
clusters in these two galaxies provide unique evidence for the view that large E galaxies can 
form by radically different evolutionary routes. 

One of the most challenging elements of IC 4051 to interpret is certainly the high 
specific frequency of its GCS. In the previous literature ( |H arris 1991| , 2000; |Blakeslee 1997 , 



1999; |Harris et al. 1998| ; [McLaughlin 1999] ) it has become conventional to associate high 



Sn with giant galaxies at the centers of rich clusters. These central BCG's or cD's can 
have had histories of star and cluster formation through inflowing gas clouds and filaments, 
mergers, and accretions (e.g., |Dubinski 1998| ) that were much more extended than for 
normal outlying ellipticals. Recently, the view has been developed that such high— Sn 
galaxies should be regarded not as "cluster-rich" but rather as "star-poor" QBlakeslee 



1997| , 1999; |Harris et al. 1998| |McLaughlin 1999]) . In this scheme, we postulate that the 



protogalactic gas started forming globular clusters at early times at a normal efficiency rate, 
but was then disrupted (perhaps by supernova-driven galactic winds, or by tidal shredding 
during infall; cf. the papers cited above) before its star formation could run to completion. 
The leftover gas now remains around these galaxies as their hot X-ray halos. This picture, 
however, assumes that the globular clusters form earlier than the bulk of the field stars 
in any given round of star formation - not an implausible requirement given the bulk of 
the observational evidence for starburst systems (see Harris 2000) and given that globular 
clusters emerge from the densest, most massive protocluster clouds. 

McLaughlin (1999) defines a globular cluster formation efficiency, measured empirically 
as the mass ratio 

M* + M gas 

where M* and M gas are the masses within the galaxy in the form of visible stars and in the 
X-ray gas respectively. He finds that e is essentially identical in the well studied Virgo and 
Fornax systems M87, NGC 4472, and NGC 1399 (despite their very different Sn), providing 
evidence for a "universal" globular cluster formation efficiency e ~ 0.26% relative to the 
initial protogalactic gas supply. The total mass ratio e is a more important indicator of 
cluster formation than Sn, which is only a measure of the cluster numbers (or equivalently 
total mass) relative to the galaxy light. In other words, Sn is a measure of only the gas 
mass M* that got converted to stars. Additional support for the near-universality of e in 
several other BCG's is found by Blakeslee (1999) . 



In this view, any high—SN galaxy should then be surrounded by a massive X-ray gaseous 
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component whether or not it is a centrally dominant giant. Notably, IC 4051 is indeed 
one of the few Coma ellipticals with an individually detected X-ray halo. Dow & White 
(1995), from ROSAT observations of the Coma core region, find that IC 4051 is detectable 
at the 2 — cr level in the soft X-ray range 0.2 — 0.4 keV, but not in the higher 0.4 — 2.4 keV 
range. If it were at the ~ 6.3 times closer distance of Virgo, IC 4051 would have a total 
Lx — 5 x 10 41 erg s _1 . This level makes it quite comparable with the Virgo giant NGC 



4472, which has Lx — 6 x 10 41 erg s 1 in the soft X-ray regime QFabbiano et al. 1992| ; jlrwin 



fc Sarazin 1996| ; |Matsumoto et al. 1997| ; |Buote fc Fabian 1998|) . However, this amount of 



X-ray gas corresponds to only ~ 5% of the stellar mass ( |McLaughlin 1999|) and NGC 4472 



as expected, has only a "normal" specific frequency level Sn — 5. 

With our adopted distance ratios for Virgo and Coma, we find that IC 4051 is about 
half as luminous as NGC 4472, so if it has a roughly similar amount of X-ray gas mass, this 
gas would only make up ~ 10% of its stellar mass. Along with Sn — 11, we find that these 
parameters convert to a present-day value for the mass ratio in IC 4051 of e ~ 0.005, twice 
as large as McLaughlin's (1999) fiducial value. 

Nominally, it therefore seems that IC 4051 acts against the paradigm of a universal 
globular cluster formation efficiency. An obvious possibility, however, is that IC 4051 
originally did possess much more gas shortly after its main era of globular cluster formation, 
but that most of this unused material was quickly stripped away as IC 4051 went through 
its first few passages of the Coma core. This gas would have joined the general reservoir of 
hot gas spread throughout the Coma potential well. The same mechanism which resulted 
in this galaxy's compact structure might then have plausibly left it with the unusual 
combination of high Sn and modest amount of X-ray gas that we now see. 

A situation which would act much more strongly to falsify McLaughlin's case for a 
universal e would be the opposite one: that is, a galaxy with a massive X-ray halo but 
a "normal" or subnormal Sn ^4. In such a case it would be much harder to avoid the 
conclusion that the formation efficiency of globular clusters was genuinely different (and 
low). Does the central Coma giant NGC 4874 present us with such a case? As we found 
in Paper II, NGC 4874 is not a high-S^ system and is embedded within a very massive 
X-ray envelope. This X-ray gas is, however, so extended that must belong to the general 
Coma potential well as a whole, with no detectable concentrated component that can be 
associated with NGC 4874 itself ( Dow fc White 1995| ). Thus there are ambiguities in the 



interpretation that are hard to circumvent. Better candidates would be E galaxies with 
massive X-ray halos that are not at the centers of rich clusters. 

Finally, we may compare the interesting case of IC 4051 with that of its Coma neighbor 
NGC 4881 ( [Baum et al. 1995| ), a giant E galaxy of similar location, size, and structure. 



- 15 - 



Curiously, NGC 4881 holds a GCS of low specific frequency (Sn 2) which appears to be 
almost entirely metal-poor, just the opposite of IC 4051. It has no significant amounts of 
X-ray gas ( pow fc White HJ§5| ). We speculate that NGC 4881 may have resulted from the 
merger of smaller galaxies in which these metal-poor globulars had already formed. These 
mergers should have been rather gas-poor to prevent the formation of newer and more 
metal-rich clusters. This is, however, an extremely sketchy interpretation, and there is an 
obvious problem with the much higher metallicity of the host galaxy light (how did the bulk 
of the giant E galaxy form at higher metallicity without leaving behind some metal-rich 
globular clusters? See Paper II for additional discussion). 

The Coma ellipticals clearly present a wide range of GCS characteristics that strongly 
challenge the array of current galaxy formation models. 
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Table 1. Completeness Function Parameters 
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Table 2. Mean Color vs. V Magnitude 
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Table 3. Mean Color vs. Radius 
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Table 4. Binned GCLF Data 
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Table 5. GCLF x 2 Fitting Results 
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Table 6. Radial Density Profile 
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Fig. 1. — Image classification plots for all detected objects on the deep V exposures. The 
radial image moment r x (see text) is plotted against V for each of the four WFPC2 CCDs. 
The dashed lines show the adopted boundaries separating star-like objects (below the line) 
from nonstellar ones (above). The vast majority of the starlike objects are globular clusters 
around IC 4051. 
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Fig. 2.— Completeness functions for the V (F606W) and I (F814W) photometry. The 
solid dots (leftmost curve) represent the PCI data, solid squares represent WF2, solid 
triangles WF3, and open circles WF4. The lines through each set of points show the Pritchet 
interpolation function curves described in the text and parametrized in Table 1. The PCI 
photometry has a noticeably brighter limiting magnitude. 
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Fig. 3.— Mean photometric uncertainty as a function of magnitude, derived from the 
artificial-star tests. The dashed line represents the PCI, while the solid lines are for the 
WF2,3,4 chips. Here and in the completeness functions, there are no significant differences 
among the three WF chips. 
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Fig. 4. — Color-magnitude distribution for all measured starlike objects with (V — I) colors. 
The vast majority of these are globular clusters in IC 4051. Objects within 80" of the center 
of IC 4051 are plotted in the left panel, and objects lying beyond 80" in the right panel. The 
solid line shows the 50% detection completeness limit in /. 
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Fig. 5. — The metallicity distribution function (MDF) for the globular cluster system in 
IC 4051. Number of objects per 0.05-mag bin, after subtraction of background (see text), 
is plotted against (V — I). The [Fe/H] scale at top follows the linear conversion relation 
given in the text, and should be taken only as schematic for [Fe/H] > 0. The best-fit single 
Gaussian function is shown, with (V — I) = 1.12 and a(V — I) = 0.13. The dashed line 
shows the color distribution for the globular clusters in the Virgo giant M87 ( |Kundu et al. 
19991) . 
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Fig. 6. — V — I colour plotted against projected galactocentric distance, for bright objects 
(V < 26.0) for which the photometry is highly complete at all radii and the measurement 
uncertainties are smallest. Horizontal lines at {V — I) = 1.17, 1.00 are drawn at the suggested 
two subpopulations; the bluer of the two components is almost absent within 10". 
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Fig. 7. — The color distribution of the clusters, fitted by a two-component Gaussian model 
as described in the text. The two subcomponents, centered at (V — I) = 1.00, 1.17 and each 
with a = 0.10, are shown as the dashed lines, and their sum as the smooth solid line. 
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Fig. 8. — Spatial distribution of the globular clusters used to define the GCLF. Dotted 
lines outline the boundaries of the four CCDs. The zone outside the outer dashed circle 
at R = 80" is used to define the "background" number density of objects. No data within 
R = 10" (inner dashed circle) were used for the GCLF analysis. Note the obvious high 
concentration of the GC population toward the center of IC 4051. 
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Fig. 9. — GCLFs (number of objects per 0.3-mag bin, after correction for incompleteness 
and subtraction of background) for the four separate WFPC2 chips. The dashed line in each 
graph represents the 50% completeness limit of the photometry, while the solid line shows 
the best-fit Gaussian function (see next Figure). 
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Fig. 10. — The GCLF for all combined V data (the sum of the four plots in the previous 
figure). The dotted and dashed vertical lines represent the 50% completeness limit of the 
PCI and WF data respectively. The solid line is the best-fit Gaussian function to the binned 
data for an assumed GCLF dispersion ay = 1.6, while the dashed line is the best-fit Gaussian 
for oy = 1.4. 
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Fig. 11.— The luminosity distribution function (LDF, or number of clusters per unit 
luminosity L/L Q ). The vertical dashed line shows the luminosity at which the GCLF (from 
the previous figure) reaches its peak or "turnover", equivalent to log(L/L©) ~ 4.8. The 
dotted line is a weighted best-fit power-law function N ~ L' a for the restricted range 
4.8 ^ log(L/L ) ^ 6.5, yielding a = 1.75. An unweighted fit to all the data (dot-dashed 
line) yields a somewhat steeper slope a ~ 2.05 affected by the downturn at the bright end. 
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Fig. 12.— Radial surface density profile for all detected objects with V < 27.0. For 
R^80" (vertical dashed line), the number density a begins flattening off to its asymptotic 
background level, which we adopt as a h = 0.02 arcsec" 2 (horizontal dashed line). 



-37- 




log R (arcsec) 



Fig. 13. — The radial projected density profile of the globular cluster system. Solid dots 
with error bars are taken from Table 7 and assume a background density = 0.02 arcsec -2 . 
The solid line is the R— band photographic surface intensity profile of the halo light, from 
Strom & Strom (1978), while the dotted line is the R— band CCD profile from Jorgensen et 
al. (1992). Both are shifted vertically arbitrarily to match to the GCS profile. The dashed 
line is the best-fit King model discussed in the text, with core radius R ~ 10" and c = 1.45. 



